New Journal of Physics

Environmental variability can strongly alter coexistence among competing species and their extinction risk, particularly when population dynamics are shaped by behavioral interactions, such as fear. In this work, we develop a novel stochastic differential equation competition model that incorporates both non-consumptive fear effects and environmental variability to investigate how behavioral interactions influence species coexistence under random fluctuations. Our result reveals that environmental stochasticity can drive species to extinction even when the corresponding deterministic system admits coexistence. In particular, under an explicit stability condition on the fear and competition parameters and sufficiently strong averaged noise intensities, we prove that both competing species become extinct exponentially almost surely. Conversely, we derive a stochastic persistence criterion in terms of fear, competition, and noise-induced suppression parameters for the fearful species. We further demonstrate that environmental noise may reverse classical competition-exclusion outcomes, leading to qualitatively different long-term dynamics from those predicted deterministically. These results provide rigorous thresholds separating stochastic extinction from persistence and highlight the critical role of environmental variability in fear-mediated competitive ecosystems. From an applied perspective, these results provide insight into how behavioral interactions and environmental variability influence species survival, with potential applications in ecological management and conservation.

Diffuse gliomas are deadly because the individual tumor cells invade - they travel far from the imageable mass, so it is impossible to remove the tumor completely. On the cellular level, glioma cells seem to be in either a "go" state (in which they do not divide) or a "grow" state (in which they do not migrate). We investigate what this tiny choice has to say about the large-scale speed of the invasion front and whether the implication is sufficiently strong to rule out the classical description of the Fisher-Kolmogorov-Petrovsky-Piskunov (Fisher-KPP) type, in which a single phenotype migrates and proliferates. We derive a two-phenotype reaction-diffusion model with density-dependent switching, and we prove the cooperative (quasi-monotone) structure and the associated comparison principle and study travelling-wave solutions of the model. A leading-edge linearization gives minimal front speed as minimizer of an explicit dispersion relation, and direct simulation verifies the predicted speed. In the experimentally relevant fast switching limit, we find a closed-form expression for the speed, that is, we obtain an effective Fisher-KPP equation with rescaled diffusivity and growth rate, with the fractions of the phenotypes. The "go-or-grow" (GoG) front can move at a maximum speed of half the Fisher speed for the same single-cell motility $D$ and proliferation rate $r$, which occurs only when the cells divide their time equally between the two phenotypes. This bound is directly testable: measurement of the front speed, plus independent determination of $D$ and $r$, discriminates the two hypotheses, and in the GoG case, yields recovery of the phenotype balance. We then extend the result to anisotropic (DTI-informed) invasion along white-matter tracts and discuss implications for understanding clinical measurements of growth rate.

Reaction diffusion (RD) systems play a fundamental role in numerous biochemical and biophysical processes. Here, we present a novel analytical framework for solving RD equations by applying the Wentzel Kramers Brillouin Jeffreys (WKBJ) formalism to Ca nanodomains generated by individual membrane channels, a widely used paradigm for intracellular Ca signaling. Previous models have primarily focused on stationary Ca nanodomains while neglecting diffusion and saturation of intracellular Ca buffers and sensors. In contrast, we derive analytical solutions without these simplifying assumptions. Our analysis demonstrates that sustained Ca influx generates continuously expanding distributions of free Ca, whereas Ca bound buffers and sensors propagate as traveling waves. These predictions are supported experimentally by measurements of one-dimensional fluorescence profiles produced by single-channel activity and two-dimensional profiles generated by whole cell Ca currents. The analytical framework developed here readily extends Michaelis Menten type kinetics to reaction diffusion systems and may therefore be broadly applicable to biochemical and biophysical processes in which diffusion cannot be neglected.

In order to feed a growing global population without silencing nature, conceiving agricultural management strategies reconciling yield and conservation goals is key. Using numerical simulations of a metacommunity model, we explore the possibilities for compromise offered by spatial management strategies of farmed areas. Each strategy is characterized by its farming intensity, the proportion of farmed lands and their spatial aggregation. We show that achieving equitable yield-biodiversity compromise is difficult. While conciliatory strategies offering top yield and biodiversity are typically not possible, accepting slightly lower yields (ie, "Pretty Good Yield strategies") allows to recover substantial biodiversity. Such reconciliation possibilities are limited for species with small dispersal. Yield increases mainly through farmland expansion, whereas farming intensity strongly influences biodiversity, increasing it at low intensity before decreasing with further intensification. Finally, we demonstrate that reconciliation is easier if agricultural production relies on biodiversity through ecosystem services.

Biomolecular condensates are dynamic assemblies whose functions depend on continuous exchange of molecular components with the surrounding environment. While scaffold molecules drive phase separation and condensate architecture, many functional components are clients that are recruited through interactions with the scaffold-rich environment. Despite their prevalence, how client-scaffold interactions shape client exchange dynamics remains poorly understood. Here, we develop a reaction-diffusion model for client exchange in scaffold-driven condensates, in which clients switch between a scaffold-bound state and an unbound state. Bound clients exchange through scaffold-mediated transport, whereas unbound clients diffuse through the pore space of the condensate. Using the fluorescence recovery of fully photobleached condensates as a measure of client exchange, we compare transport through these two pathways with bound-unbound conversion and identify three limiting regimes. In the slow-conversion regime, bound and unbound clients recover through distinct scaffold- and pore-mediated pathways. In the intermediate-conversion regime, recovery of bound clients becomes limited by client unbinding. In the fast-conversion regime, local equilibrium between bound and unbound clients produces an effective single-state recovery. We further propose a unifying description that connects these regimes and quantitatively captures the apparent recovery timescales extracted from numerical simulations across condensate sizes. Our results provide a framework for interpreting component-specific exchange dynamics, and highlight client size, client-scaffold binding, and condensate porosity as key regulators of client turnover in multicomponent condensates.

Directed migration along chemical gradients controls immune surveillance, development, and cancer invasion. However, the same chemical cue can produce different responses depending on its concentration and whether cells move alone or in groups. For example, in steep gradients, isolated malignant lymphocyte cells migrate away from the chemoattractant source, whereas clusters of the same cells continue to migrate toward it. Here, combining computational modeling and experimental observations, we show that this reversal is governed by coupled mechanisms acting across molecular, cellular, and collective scales. At the single-cell level, our model predicts that receptor endocytosis generates a feedback that produces a nonmonotonic surface receptor density with increasing chemoattractant concentration. Above a critical concentration that depends on the cell's volume-to-sensing-area ratio, receptor depletion reverses cell polarity and drives chemorepulsion. However, in clusters, cell-cell contacts reduce the membrane area exposed to ligand, increasing the volume-to-sensing-area ratio, thus increasing the critical concentration and preserving chemotaxis. An agent-based model incorporating these mechanisms quantitatively reproduces the sign reversal of the migration index across gradient steepness and cluster size. We show that collective rearrangements further stabilize chemoattraction with exchanges between the cluster rim and core helping remove chemorepulsive cells from the leading edge, keeping their fraction below the threshold required to reverse cluster migration. The model further predicts, and experiments confirm, that increasing ambient ligand concentration while keeping the gradient fixed reduces cluster chemoattraction. Our results identify receptor trafficking, cell geometry, and cluster fluidity as physical determinants of collective directional decision-making, with implications for immune cell homing, tissue morphogenesis, and cancer dissemination.

Cell shape and mechanics are intricately connected and tightly regulated by mechanochemical events including biochemical signaling, cytoskeletal remodeling, and plasma membrane mechanics. While experimental advances in microscopy have shed light on the intricate coordination involved in cell shape change in response to different cues, the ability to conduct three-dimensional simulations in realistic geometries remains an open computational challenge. In this work, we develop a finite-element framework that incorporates advection-diffusion-reaction equations coupled with equations governing the kinematics of a deformable interface representing the cell membrane. We applied this framework to three distinct coupled mechanochemical systems, each governed by geometric partial differential equations, resulting in large deformations of the interface. In all three examples, our simulations revealed the emergence of feedback between cellular signaling, cytoskeletal organization, and cell shape. In our first two sets of simulations, we observed that cell migration and neutrophil protrusion were regulated by membrane tension-mediated feedback. In our final application, we predicted shape changes of a dendritic spine starting from a realistic geometry, and found that the complex shape of the spine gives rise to localized regimes of actin cytoskeleton remodeling not previously observed with idealized geometries. Thus, our finite-element framework allows us to generate new mechanistic insights for biophysical problems.

Under complex perturbation regimes, biodiversity dynamics show temporal variability in species and community abundance around long-term population trends. Many species indeed show long-term declines while other species increase, putting natural communities far from stationary regimes, while variability is often studied near equilibrium. We contribute to bridging this gap by investigating population and community variability during long-term trends caused by press perturbations in stochastic models of population dynamics. By estimating the deterministic changes in mean and variance during the transient regime, we show that population variability deviates from stationary expectations. Moreover, the deviation strongly depends on the sign of the population trends: increases generate excesses of variability while declines generate deficits. Scaling up to community variability, we propose a decomposition of community variability deviation, allowing to highlight that community variability in the transient regime depends on how the press perturbation is distributed within species relative abundances and growth rates. These results challenge the equilibrium assumption and open new perspectives for the study of the variability of ecological systems under multiple perturbation types.

Auxetic metamaterials exhibit negative Poisson's ratio behaviors due to their architecture of periodically arranged unit cells. Although mechanical metamaterials are well established at the macroscale, programmable auxetic units remain scarce at the nanoscale. DNA origami offers a promising platform to bridge this gap, but design principles for dynamically deformable 3D auxetic nanostructures remain largely unexplored. Here, we develop design strategies for such 3D auxetic metastructures built from wireframe DNA origami. As a model system, we use a 3D re-entrant triangular unit composed of double-stranded DNA (dsDNA) bundle edges connected by single-stranded DNA (ssDNA) joints. Using coarse-grained molecular dynamics (MD) and umbrella-sampling free-energy simulations, we examine how edge design and joint-connection scheme govern auxetic responses and the energetics of the structural transformation. Our results show that auxetic performance and deformation energetics emerge from the coupled effects of DNA bundle rigidity and connector mechanics at the joints. This study provides mechanistic insights and design guidelines for programmable auxetic motion and energetics in 3D DNA origami metamaterials, advancing the development of stimuli-responsive nanomechanical devices.

There is a vast literature on the voltage gating of ion channels, with a fairly large fraction concerned with potassium channels, especially of the KV1 family, including Shaker. Experimental evidence derived from protein structure has been interpreted to give gating mechanisms that largely disregard water. We propose that the K+ ion, in order to pass through the gating region and enter the cavity pore, must be largely dehydrated. Competitive interactions of each single hydration shell water at the gate, with counterions, protein, or other water molecules, can remove one water at a time. There are several such interactions for the ion hydration shell; for the ion to pass through the gating region, there must be enough such interactions to leave the ion with at most two hydrating water molecules, in which case the gate is open. Protein conformational changes are secondary, small, and mostly unimportant. The hypothesis has a second part: protons, previously shown to be candidate carriers of the gating current (Kariev and Green, JPC B, 2019, Membranes, 2022, 2024) are capable of reaching the gate; adding four protons to the gate prevents dehydration, leaving the ion with at least three hydrating water molecules, enough to block passage. Quantum calculations presented here support the dehydration part of the hypothesis; they also mostly support the second part, concerning the protons, but further work will be required to fully confirm this. The hypothesis explains the experimental finding that the P475D mutant is essentially constitutively open, while the P475S mutant, with a wider gate opening, is closed at all relevant potentials; the computations presented here show the mechanism for this in detail, further confirming the first part of the hypothesis, and largely but not completely confirming the second part, concerning protons, while showing where further work is needed. This mechanism can also qualitatively account for flicker noise and fluctuations, and their consequences.

Every brain must adapt to an unpredictable world, yet individuals differ in how readily they learn. Theoretical work suggests that learning is fastest when a system, whether biological or synthetic, is initialized in a state close to instability - i.e., near criticality - because critical dynamics are imbued with a diverse repertoire of patterns and multi-scale correlations. Here, we empirically estimate distance to criticality in the brain and show that it predicts the rate of adaptability underlying learning, neuronal tuning, and general intelligence. In mouse motor cortex, proximity to criticality forecasts learning rate of two future complex tasks: prey capture hunt and ladder crossing. In contrast, distance to criticality predicted neither an animal's naive ability nor its asymptotic skill - isolating the rate of learning itself. In visual cortex of young ferrets, proximity to criticality predicts how strongly experience reshapes neural tuning. In human frontal cortex, it correlates with general cognitive ability. A minimal recurrent network model reproduced these results and offers a mechanism: proximity to criticality defines the timescale over which a system can learn from its past experiences, directly setting the rate of learning. A single dynamical property can account for the capacity to learn, from artificial networks to the mammalian brain.

Single-nucleus transcriptomic atlases offer an unprecedented opportunity to connect cellular molecular states with Alzheimer's disease (AD) neuropathology, but whether these profiles encode reproducible, predictive information about pathological burden remains unclear. We present the SEA-AD DREAM Challenge, an open, international, model-to-data competition built on the Seattle Alzheimer's Disease Brain Cell Atlas to predict Alzheimer's disease neuropathological severity from single-nucleus RNA-sequencing data. Participants developed containerized models to predict categorical neuropathological staging, including overall Alzheimer's disease neuropathologic change, Braak stage, Thal phase, and CERAD score, as well as quantitative amyloid-{beta} and phospho-tau burden measured by 6E10 and AT8 immunohistochemistry. Across 17 eligible teams from 15 countries, the crowdsourcing framework enabled systematic comparison of diverse computational approaches and surfaced a broad landscape of modeling strategies and candidate predictive features. Top-performing methods achieved near-perfect prediction of categorical staging, with the best submission reaching a quadratic weighted kappa of 1.0 for the Overall AD Neuropathological Change score (ADNC), and competitive prediction of quantitative pathological burden in held-out data, with a best concordance correlation coefficient of 0.48. Post hoc perturbation analyses revealed that top categorical-stage predictions relied heavily on donor-level metadata-driven signals rather than transcriptomic features, whereas quantitative pathology prediction was more robust and supported by transcriptomic and cell-type-associated features with potential biological relevance to AD progression. The challenge also introduced the first AI Agent Track in a DREAM Challenge, providing an early benchmark for autonomous and human-guided agentic model development in single-cell neuroscience. This work demonstrates that single-nucleus transcriptomes encode substantial information about Alzheimer's disease pathology, establishes a reproducible benchmark for molecular neuropathology prediction, and highlights critical principles for designing privacy-preserving, leakage-aware community challenges using deeply phenotyped human brain data.

Temporal niche partitioning is a strategy for reducing interspecies competition and strengthening reproductive isolation. It relies on animals confining their daily activity to distinct diurnal, crepuscular, or nocturnal windows. However, a hardwired temporal niche is only advantageous under stable, predictable ecological regimes; surviving dynamic environments demands behavioral flexibility. Yet, it remains unclear how animals override rigid biological constraints to rapidly exploit transiently available fitness-critical time windows. To address this, we leveraged the twilight-active, species-rich Drosophila genus and monitored their daily activity under naturalistic conditions. Here, we show that intense sociosexual interactions rapidly drive a species-specific reformatting of their canonical crepuscular niche. The dominant sensory modality used for sexual communication predicts niche shift direction: reliance on chemosensation for courtship redirects behavioral activity into the night, while visual reliance shifts it into the day. This temporal plasticity bypasses the circadian clock, instead operating via a conserved dopaminergic pathway. Dopamine operates a dual-output brain circuit that simultaneously inhibits sleep and sustains sexual motivation. Our results reveal how mating imperatives decouple behavioral timing from circadian command, enabling conditional colonization of otherwise restricted temporal windows. Ultimately, by driving the divergence of previously overlapping niches, sociosexually induced temporal plasticity provides a powerful mechanism for sympatric coexistence in crowded environments.

Extrahepatic delivery of lipid nanoparticles (LNPs) to non-phagocytic cells is a major challenge, with the leading strategy involving surface functionalization with target-specific monoclonal antibody (mAb) ligands. We investigate the stability of mAb-conjugated LNPs using two anchoring systems: the commonly used DSPE-PEG2kDa-maleimide and a block copolymer, PCL5kDa-b-PEG2kDa -maleimide, with the hypothesis that conjugation to a 150,000 Da antibody could overwhelm the relatively small ~600 Da aliphatic anchor on the PEG-lipid in vivo. Shedding of the mAB would compromise targeting. Conjugation integrity following IV injection was assessed by tagging LNPs and mAbs with metal ion tracers that could be quantified by ICP-MS. Results show that DSPE-PEG-mAb rapidly (within 1h) dissociates from LNPs in blood, leading to accelerated LNP clearance. In contrast, mAbs conjugated using PCL-b-PEG remained stably associated with the LNP over the 24h circulation and clearance of the construct. Results are connected to a thermodynamic model that reproduces experimental findings for PEG-anchor(-mAb) shedding in vitro and in vivo. This study identifies anchoring strength as a critical, unconsidered parameter for in vivo performance when conjugating mAbs to LNPs for extrahepatic delivery.

The incidence of early-onset colorectal cancer (CRC) has risen sharply in recent decades1, yet the biological basis underlying the distinct behavior of tumors arising in young versus aged tissues remains poorly understood. Here we show that aging reprograms the epigenetic landscape of the colon, restricting colon tumor growth through stable silencing of developmental and fetal gene programs. We find that colon tumors arising in aged mice are intrinsically less proliferative than those arising in young animals. Multi-omic profiling of normal colon and colon tumors reveals that aging drives DNA hypermethylation, loss of Polycomb-associated chromatin states, and reduced chromatin accessibility at a defined set of developmental genes that are bivalent (marked by both H3K27me3 and H3K4 methylation), transcriptionally active in colon tumors from young animals and repressed in both tumors and normal tissue from old animals. Among the genes most strongly repressed in old animals is Tacstd2 (Trop2), a regulator of fetal intestinal programs and epithelial stemness. Pharmacologic inhibition of DNA methylation reactivates the aging-silenced gene network in organoids from old animals, whereas genetic disruption of Tacstd2 suppresses growth and developmental transcriptional programs in young tumor organoids. TACSTD2, fetal gene signatures, and the aging-associated bivalent gene program are likewise repressed in late-onset vs. early-onset human colorectal cancers. Collectively, these findings identify age-associated epigenetic silencing of developmental gene programs as a causal mechanism that constrains colorectal tumor growth and provide a mechanistic framework for understanding the distinct biology of early-onset colorectal cancer.

The spatial arrangement of immune cells in the tumor microenvironment (TME) varies widely, from dispersed to clustered and tumor excluded to infiltrating. Multiplexed spatial profiling is an effective means of characterizing tumor-infiltrating lymphocytes (TILs) and immune complexes such as tertiary lymphoid structures (TLS) in the TME. However, few approaches have been described for objectively parametrizing patterns of immune organization and assessing their association with biological or clinical variables. This makes it difficult to evaluate whether a set of tumors is relatively immunologically cold or hot. Here we describe an intuitive set of statistical tools (available in the R package, tlsR) for characterizing lymphocyte patterns in the TME of solid cancers. We apply tlsR to primary prostate cancer (PCa), which is often described as immunologically cold. Using a cohort of 29 radical prostatectomy specimens stratified into low Gleason-grade (LGG; n=15) and high Gleason-grades (HGG; n =14) we show that HGG PCa is significantly more infiltrated than LGG PCa with lymphocytes organized into B cell or T cell enriched immune clusters (BICs and TICs). A subset of these ICs have the B and T cell zonation and follicular dendritic cells characteristic of a bona fide TLS. HGGs are also enriched with ICs containing precursor exhausted T cells (Tpex) and proliferating B cells and their tumor compartments harbor granzyme-B+ cytotoxic T cells in contact with cancer cells. Thus, far from being cold, a subset of HGG PCa has features associated with active immune surveillance, a finding with implications for emerging PCa immunotherapies.

Blockade by sedimentary particles, such as mineral crystals, is a continuous risk the kidney tubule faces. To prevent that, kidney resident macrophages form transepithelial protrusions and remove intratubular sedimentary particles, a behavior particularly prevailing in the medulla over the cortex. However, the molecular mechanisms underlying this characteristic behavior of medulla macrophages are incompletely understood. In this study, we identified that the medulla had higher mechanical stiffness than the cortex in steady state, which was further elevated when kidney stone formed. Increased tissue rigidity was sensed by medulla macrophages via mechanoreceptor Piezo1, which promoted macrophage protrusion formation and their ability to clean the tubules. Loss of Piezo1 expression in kidney macrophages predisposed mice to intratubular accumulation of mineral crystal in steady state and accelerated kidney stone formation during oxalate intake challenge. Signaling via Piezo1 mobilized molecules involved in cell adhesion and protrusion assembly, including Talin2 and focal adhesion kinase (FAK). Finally, we developed a first-of-its-kind cell-based therapy for the treatment of experimental nephrolithiasis by exploiting macrophage Piezo1 activity, and this strategy shows great promise for future translational research.

Spectral bleed-through remains a persistent practical problem in multichannel fluorescence microscopy. Signal from one fluorophore can be recorded in the detection channel of another, thereby biasing intensity measurements, inflating apparent colocalization, and complicating the interpretation of dynamic microscopy data. Although many correction strategies exist, routine workflows often remain fragmented across ad hoc scripts, manually tuned graphical procedures, or method-specific blind-unmixing implementations with limited provenance. Here we present spectral-unmixing, an open-source Python package for reproducible linear spectral unmixing in multidimensional microscopy stacks. The package unifies directed two-channel correction with multiple alpha-estimation strategies, optional bidirectional two-channel correction through explicit inversion of a 2 x 2 mixing model, and PICASSO-family blind unmixing for multichannel data. Microscopy inputs are normalized at the API boundary to canonical TZCY X stacks, allowing the same unmixing code to be applied across file formats without manual axis handling. Machine-readable sidecar reports preserve the effective processing configuration and estimated coefficients for every output, so that workflows can be audited and reproduced. Synthetic and real-data-derived benchmarks show that the implemented workflows accurately estimate and correct bleed-through when their model assumptions are satisfied. In fixed-alpha two-channel simulations, the mean-ratio and linear-fit estimators recovered {approx} 0.283 for a ground-truth value of 0.28 and reduced target-channel normalized root mean squared error from approximately 0.029 to 0.003. In time-varying simulations, per-time-point estimation tracked coefficient drift substantially better than reference-time-point estimation. Bidirectional inversion recovered reciprocally mixed channels accurately when coefficients were known or well estimated. PICASSO-family benchmarks further showed a practical trade-off between reducing residual inter-channel dependence and preserving fluorophore identity, with MATLAB-style workflows behaving more conservatively and source-sink formulations providing stronger dependence suppression when meaningful directional priors are available. Together, these elements make spectral-unmixing a practical, transparent, and extensible platform for reproducible spectral unmixing of fluorescence microscopy data in neuroscience and other quantitative bioimage-analysis settings.

We showed before that the endogenous peptide Nociceptin blocks the premature differentiation of oligodendrocytes (OLGs), preventing untimely precocious myelination in the developing brain. Consistent with this early function, Nociceptin brain expression is developmentally regulated, sharply decreasing with the initiation and progression of myelination. However, we now found that at difference with controls and relapsing-remitting multiple sclerosis (RRMS), Nociceptin levels are highly elevated in cerebrospinal fluid from patients with the most severe progressive MS (PMS) forms. This questioned whether Nociceptin early developmental effects could be latter recapitulated, interfering with remyelination in PMS. This possibility was tested by inducing experimental autoimmune encephalomyelitis in older mice, at an age equivalent to that with increased risk of RRMS transition into PMS. Older animals develop persistently highly debilitating clinical symptoms, and display both brain and spinal cord demyelination. Importantly, these mice exhibit elevated brain Nociceptin levels, and their treatment with an antagonist of the Nociceptin receptor (NOR) elicits a regression of clinical scoring that is accompanied by higher ratios of OLGs/OLG progenitor cells, increased myelination, and reduction of reactive astrocytes. These findings suggest that Nociceptin may be a crucial player in the age-related progression of MS; interfering with OLG maturation and remyelination, and perhaps further exacerbating neurological dysfunction by targeting astrocyte populations. The upregulation of Nociceptin secretion by human astrocytes in response to proinflammatory cytokines, also points to this peptide as a mediator of microglia-astrocyte interactions supporting MS progression with aging. NOR may offer a novel pharmacological target for ameliorating the devastating effects of MS progression.

Recent advances in spatially resolved transcriptomics (SRT) enabled measurement of sets of pathway genes activity within tissues. However, existing gene set activity scoring methods overlook spatial dependencies among tissue locations, restricting their ability to capture region-specific pathway activities associated with disease pathology or cellular communication. Moreover, these methods lack significance-level inference for activity scores, provide limited interpretability of gene-level contribution to a pathway, and scale poorly to advanced large-size SRT datasets. To address these limitations, we present GESSO (Gene sEt activity Score analysis with Spatial lOcation), a spatially informed gene set scoring method adaptable to diverse SRT platforms. GESSO models gene set activity levels through a graph-regularized matrix decomposition algorithm, jointly inferring spatially coherent gene set activity scores (GASs) and interpretable metagene weights that capture gene-level contributions. It further implements a permutation-based local significance test and a stratified low-resolution approximation that scales to high-resolution SRT datasets such as Visium HD, Stereo-seq, and Xenium Prime. Across 13 datasets from five SRT platforms, GESSO outperformed all existing methods in accuracy, calibration, interpretability, and scalability. Applications revealed novel biological programs, including spatially confined EMT activation within tumor-stroma interfaces, developmental signaling gradients across embryonic tissues, and coordinated B-cell, T-cell, and signaling pathways within germinal centers of human lymph node tissue, revealing the spatial organization of immune function at subregional resolution.